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ABSTRACT 

Viscous resistance to changes in the volume of a gas arises when different 
degrees of freedom have different relaxation times. Collisions tend to oppose 
the resulting departures from equilibrium and, in so doing, generate entropy. 
Even for a classical gas of hard spheres, when the mean free paths or mean 
flight times of constituent particles are long, we find a nonvanishing bulk 
viscosity. Here we apply a method recently used to uncover this result for a 
classical rarefied gas to radiative transfer theory and derive an expression for 
the radiative stress tensor for a gray medium with absorption and Thomson 
scattering. We determine the transport coefficients through the calculation 
of the comoving entropy generation. When scattering dominates absorption, 
the bulk viscosity becomes much larger than either the shear viscosity or the 
thermal conductivity. 



1 INTRODUCTION 

Both the thermal and mechanical effects of radiation on matter play a significant role in the 
fluid dynamics of hot objects, from the early universe to early stars. The thermal effects may 
be reasonably well accounted for by the use of the Eddington approximation (|Unno fc Spiegel 
1966|) , though special circumstances may call for improved descriptions. In particular, when 
travel times of photons are not negligible, retardation effects may be significant ( pelache~fc 
Froeschle 1972|) , as indeed they are when we go to treat the dynamical effects of the radiation 
field. 

An early description of the dynamical effects of radiation on matter was provided by L. 

© 2000 RAS 



2 X. Chen and E. A. Spiegel 

H. Thomas ( Thomas 1930|) who expanded the specific intensity in terms of the mean free 
path of photons, as Hilbert had previously done for the one-particle distribution in kinetic 
theory ([Kogan 1969|) . Each of them assumed that the relevant mean free path is much smaller 
than the characteristic macroscopic length scales of the material medium. Though Thomas's 
treatment included only emission and absorption, Masaki later extended it to include the 
effects of Thomson scattering ( |Masaki 1971| ; [Hsieh fc Spiegel 1976| ) and Compton scattering 
( IMasaki 19811 ). 



The effectiveness of the dissipation in those studies was characterized by the radiative 
transport coefficients — shear viscosity and conductivity — about whose evaluation there 
has been general agreement since Thomas's derivation of the stress tensor for a matter- 
radiation mixture with short photon mean free paths ( |Weinberg 1972t |Mihalas fc Mihalas 
1984). But in recent years, particular attention has been given to the bulk viscous effects, 



which are of interest for applications to cosmology and, for that case, different expressions 
have been proposed. 

A simple relativistic gas has bulk viscosity flStewart 1972|) because the ultrarelativistic 
particles and particles with small energies respond to volume changes at different rates 
and so are driven out of equilibrium with each other under expansion. Analogously, in 
radiative fluid dynamics, the temperatures of the matter and the radiation may not be 
equal in nonequilibrium processes. This too can be a cause of a bulk viscosity ([Weinberg 
1971; Mihalas fc Mihalas 1984), even though photons are ultrarelativistic. Weinberg gave 



an expression for the radiative bulk viscosity in the Thomas case with only emission and 
absorption ( |Wcinbcrg 1971| ). Subsequent work ( [Anderson fc Kox 1977] ; [Stewart 1972] ), also 
for media with only emission and absorption but based on other methods, gave different 
expressions for the forms of the transport coefficients. However, those results have been 
shown to be equivalent to Weinberg's in the first order in photon mean free path. 

Additional studies have included scattering but, in those more complicated calculations, 
the bulk viscosity has either not been derived or not extracted explicitly ( |Masaki 1971 



Hsieh fc Spiegel 1976| ; |Straumann 1976| ; [Masaki 1981| ; |Thorne 198 1|) . Hence our aim here 
is to determine the bulk viscosity coefficient from a calculation of the entropy generation 
in the comoving frame of the matter, once we have obtained suitable approximations for 
the radiative stress tensor. Since the previous studies have been limited to cases where the 
mean free paths of photons are much smaller than the prevalent characteristic scales, the 
problem of treating situations where transparent regions form in the medium has been left 
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open. This issue must be confronted at the edges stars and disks, in intense turbulence and 
in photon bubbles. A case of particular interest arises in models of the early universe when 
the mean flight times of photons are comparable to the age of the universe. 

The difficulties posed by long mean free paths are less severe for thermal than for dynam- 
ical problems since, as we mentioned, the radiative smoothing of temperature fluctuations 
is reasonably well described by the Eddington approximation, both in optically thick and 
thin regions. A reason for this is that the Eddington approximation represents a summation 
of contributions from terms of all orders in a mean free path expansion ( |Unno fc Spiegel 
|1966| ). However, since the Eddington approximation leads to a diagonal stress tensor, it fails 
completely to describe shear viscosity ( |Anderson fc Spiegel 1972] ) . Moreover, it does not give 
a good representation of the bulk viscosity that opposes changes in volume. 

Methods based on truncation of the moment hierarchy at higher moments than the 
pressure tensor do lead to shear viscosity flThorne 1981| ; (Btruchtrup 19971) , but they give 
complicated equations and are difficult to use. Another approach is to use a more general 
resummation procedure than that leading to the Eddington approximation, one that provides 
a nonlocal viscous shear tensor. Though this has been done for the case of pure absorption 
and emission (|Chen fc Spiegel 2000|) , it has not been carried out with scattering included. 
Nor is this issue to be considered here. Our purpose in this paper is rather to confront the 
problem of bulk viscosity for a medium where the photon mean free path is not short and 
where we also have to deal with scattering processes. For this, we shall adopt a procedure 
that has recently been introduced in classical kinetic theory ( |Chen et al. 2000| ) to generalize 
the usual Navier-Stokes equations. The generalized equations exhibit resistance to volume 
changes when the particle mean free paths are long, even for a gas of hard spheres. 

To test the results obtained with the new method, we have calculated the thicknesses of 
shock waves and the propagation speeds of ultrasonic waves for classical gases. The results 
found are in good agreement with experiment, whereas those obtained from the Navier- 
Stokes equations do not fare well. Though the limits of validity of the new approach are 
still being probed, the results seem good enough to make the application to the problems of 
radiative fluid dynamics seem worthwhile. Since the methods that have been used heretofore 
to compute the radiative dissipation terms parallel the Chapman- Enskog procedure flKogan 
1969| ; |Cercignani 1988|) , we may expect to encounter differences from previous results in this 



case as well. 

In what follows, we describe an expansion in photon mean free path for solutions of the 
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transfer equation, including the effects of Thomson scattering, and develop a stress tensor 
for use in media whose photon mean free paths need not be small. From this, we compute 
the rate of comoving entropy generation and provide a formula for bulk viscosity. In a later 
paper, we shall use this result to estimate the entropy generation in the standard model of 
cosmology. 



2 TRANSFER THEORY 
2.1 The radiation field 

We consider radiative transfer in a fluid medium with density p and velocity u^, each de- 
pending on location in space-time, x M , where fi — 0,1, 2, 3. There are two reference frames 
of basic interest in this work, an inertial frame, or system frame, and the comoving frame of 
the matter. In the system frame, we have 

^ = 7 (l,v), 7 = (l-v)- 1 / 2 (1) 

while in the comoving frame = (1, 0). Here, as in the following, we assume units in which 
the speed of light in vacuo and Planck's constant are each unity. We also adopt the signature 
in which the Minkowski metric is diag(l, —1,-1,-1). 

The four-momentum, p^, of a photon satisfies p^ = 0, so that we may define a basic 
null vector, 

r^=?/>=(l,n), (2) 

that characterizes the direction of the photon's motion in spacetime. With v as the frequency 
in the inertial frame, the rest frequency (as seen comoving locally with the medium) is given 
by the Doppler formula 

v = u^p^ = 7^(1 — v ■ n) . (3) 

The coordinates of the phase space of photons are x M and which we shall sometimes 
denote as x and p for brevity. We shall describe the density of photons in phase space by the 
one-particle distribution function, f(p,x). This quantity is a scalar and its integral over the 
invariant volume in momentum space gives the photon number density in space. To perform 
this integration, we introduce the invariant volume in phase space, dP = ududQ, where dQ 
is the element of solid angle ([Landau fc Lifshitz 1984j ). 
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Radiative Bulk Viscosity 5 

A principal quantity of interest in this work is the stress-energy tensor of the radiation 
field, 

= J p^p u f alP = J J n"n u l dQalu , (4) 

where the last term on the right introduces the specific intensity, defined as 

l=v 3 f. (5) 

A revealing way of writing the stress tensor in terms of basic moments of the radiation 
field is obtained by defining the directional vector with respect to the moving matter, 

p = - , (6) 

so that 

l^F = -1 and = • ( 7 ) 

Hence Z M = (0, 1), where 1 is a unit three- vector in the spatial direction of photon motion in 
the comoving frame. The angular moments of the radiation field are then 

£= fldfl , J* = [ Xl^dVt , = f Xin u dVt , (8) 



and their frequency-integrated forms are 



E= Sdu F^= T^dv = / V^dv . (9) 

Jo Jo Jo 

These moments together make up the energy-momentum stress tensor which may now be 
written as 

= Eu»u u + F»u v + F u u fJ ' + . (10) 
2.2 The Transfer equation 

When both absorption and scattering occur, the transfer equation takes the general form 
mon 19631) 



jfd.f = p(a -(3f)+pJ M(p,p')f(p')dP' -pj K(p\p)f(p)dP' (11) 

where a and (3 are scalars characterizing absorption and emission and the kernel 3?(p',p) 
is the differential cross-section for scattering a photon from four-momentum p' into p. For 
brevity, we have not indicated the dependence of / on x in the transfer equation and have 
omitted the effects of gravity. We here consider only the case of Thomson scattering for 
which we have ( |Masaki 1971| ; [lucker 1975| ) 
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U(p,p') = ? [1 + (1-1') 2 ] 8{v-i?)^ (12) 

where a is the Thomson cross-section. 

We may factor / out of the second integral of (|TTD and note that with i = 1, 2, 3 we have 
r ■ ■ An ■ ■ 

jrPdtt = -—5 tJ (13) 

where 5^ is the Kronecker symbol. Then the last term in (|IT ) reduces to —per/. Since / 
is a scalar, we may use (|5|) to find the transformation law for the specific intensity. Also, 
we may compare the covariant form of the transfer equation with the usual one for the 
comoving frame, as Thomas (Thomas 1930|) did, to find (3 = uk(u), where k is the absorption 
coefficient, and a = v~ 2 j{v\ where j is the emissivity. Since a, (3 and v are scalars, we then 
know the transformation rules for j and k once we introduce the (relativistic) Doppler 
formula (§). 

The equation of transfer may now be written as ([Hsieh fc Spiegel 1976| ) 
3(3" 

en p I p = kS + l p LV pa ) , (14) 

l07T 

where 

S = v 3 a/(3 , (15) 

e=[p( K + ( r)]~ 1 , (16) 

k=—, a = ^- (17) 

K + a K + cr 

and a comma with subscript p indicates differentiation with respect to x^. We have replaced 
SijV^ by S pa V pa , to which it is equal, for cosmetic reasons. 

Expressions for the moments that appear in the transfer equation may be computed 
from the equation itself. In doing this, we shall assume that the medium is grey, although 
we could save appearances by introducing suitable mean absorption coefficients. We shall 
not follow that practice here since it leads to serious complications in moving media, such 
as tensorial absorption coefficients, and has not repaid the effort involved. 

Next we derive expressions for the monochromatic moments defined in (Rf) by taking 



suitable moments of (H). For this purpose, we define 

h pa = -— f l p l a dQ , (18) 



47T 

which, in terms of the basic tensors of the problem, may be expressed as 
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Similarly, we have 

J dn = ^r(h^h pa + w p k + Kh v p ) . (20) 

Then, when we integrate ([[4]) over all solid angle and note that h pu V pu = —£, we find 
that 

£ = AnS-^r f n p l p dtt . (21) 
Further, on applying the second definition in (||) to (|i4|), we get 

= -e [ l p n p X )P dVt . (22) 



Finally, on multiplying flTjj ) by l p l u and integrating over angle we find that 

v „u = _ 407TK _ _ _10£^ f 

3(10 -a) 10 -a 10 - a 7 ,p v 7 

With the moments expressed in this way, we obtain the expression 

S + LlaV pa = —S - 2£ } 5 + a \ [ n p l p dn - -^rU„ / l p l v n p X p d£l . (24) 
p 3 k(10-o-)7 ,p 10-ct p 7 ' p V ; 

Hence we may write the transfer equation in the compact form 

1 = S — eCl (25) 
where we have introduced the linear operator 

C = n p d p + Cl JdVL n p d p + c 2 l p l v J dSl Fl u n p d p (26) 
with d p = d/dx p , 

Cl = J^±^L (27) 

87r£(10-a) V ' 

and 

C2 = Mlo^I) ■ < 28 > 

When there is no scattering, <t = and the foregoing equations reduce to those for the case 
of pure absorption. 



3 THE RADIATIVE PRESSURE TENSOR 
3.1 The Expansion Procedure 

We expand the intensity I as 

oo 

J = ]T I (m) e m . (29) 

m=0 
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From X( m ) , we can calculate the coefficients of the expansions of the monochromatic moments 
and their frequency-integrated counterparts from the definitions 

£(m) = J Z(m)dn , T{ m) = J iwFdn , vfo = J i (m) n u dn (30) 

and 

poo poo POO 

E(m) = I £(m)dv Ff m) = T[ m) dv Pfc = V^du . (31) 

J (J J (j J (j 

When we introduce the expansion fl2"9|) into the transfer equation fl25|) and demand that 
the expanded equation is satisfied term by term, we get 

J (0) =5, (32) 

as in comoving local thermodynamic equilibrium. For the higher orders with m > 1, we 
obtain 

1(m) = — £Z(m-l) (33) 

where C is defined in (|26|) . 

The problem studied is really one of a mixture of radiation and matter and a more 
complete description would begin with coupled transport equations. In the present version, 
we simply proceed as if the material properties are given, as did Thomas ( [Thomas 1930| ). 



In doing this, we consider that the properties of the medium are expressed, not as functions 
of x p , but in terms of the basic fields, T, 9, n p . Therefore, in the simplest case, where we 
operate on a function of T, v and n p only, we may write 

d„ = T tli d T + i> tli d s + n^d n « . (34) 

On recalling that v = and that p p does not depend on coordinate, we obtain 

v tfX = p u u^ . (35) 

Similarly, on differentiating p p = n p v and making use of ( |35l) and the identity u^u^ p = 0, 
we get 

n p y = - n p n p u P)U . (36) 
So ( |34"D takes the form 

<9 M = T^d T + i>n p u Ptfl du - n a n p u p ^ d n a . (37) 

If we further assume that, like the Planck function, S depends on v and T only through the 
combination u/T, we deduce from ( |37[ ) that 
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S :P = (T p -TnP Up ^)d T S. (38) 

We note that in the higher orders some of the derived quantities will depend explicitly 
on the derivatives of T and u p as well as on the fields themselves. Thus, more generally, 
when we go to higher order developments, as we did in computing the shear viscosity ( |(Jhen| 
|fc Spiegel 2000| ), we must allow for a dependence on those derivatives. However, here we go 



only to first order in photon mean free path and that complication does not arise; ([J7|) is 
therefore sufficient for our purposes. 

Returning to the calculation of S, we see from fl2lD that £(o) = 47r<S and that, for m > 1, 



K 



I f n p d p l {m „ l)( m . (39) 



Also, from (P2"D we have = and, for m > 1, 



Hn) = ~ / ^ ^(m-1) dn . (40) 



(m) 

Then from (g3j) we find Pgj = -^-Sh^ and, for m > 1, 

% = "I^T^M^ " 77f^ / W dn ■ ( 41 ) 

3.2 The First-Order Stress Tensor 

From the recursion relation ([33|) we have 

Z (1) = £Z (0 ) = £5 . (42) 
Then from ([3~8"D and (|39|) we get 

^ (1) = --(^ + T)5 T (43) 
where 

= (44) 
and T = u^T^. On integrating over frequency, we find 

Alt 4 

E (1) = —(-S9 + S), (45) 
where 

POO 

S= / Salu = aT i (46) 

JO 

and 5" = u^S tfl . 

In a similar way, we find for the frequency-integrated first-order radiative flux, 
F^ = ^[S >p -ASu p ]h^, (47) 
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while for the frequency-integrated pressure tensor in first order, we obtain 
40tt 

/>'" /-.. , . 

3 



^) ~ 10 - a 



"\a h^ v .AS ^ 



40vr 



where 



.pupa 



We may then introduce (15) and so write the pressure tensor to this order as 



An 

pV>v — 

3 



S-~(~S6 + S) 

K 3 



?pV 



where 



^pv 



jlT 



pupa 



U 



p.cr 



(48) 



(49) 



(50) 



(51) 



is the shear tensor and the shear viscosity coefficient is that found by Thomas and Masaki, 
_ 32nS 

* ~ 3p(Wn + 9a) ■ (52) 
As mentioned, we have elsewhere computed a nonlocal expression for the shear tensor for 
the case of pure absorption ( Chen fc Spiegel 2000|) and expect that a similar generalization 
may apply in the case with scattering, though we have not yet carried out the necessary 
calculations to confirm this. 

From (f47|), we may identify the thermal conductivity as the coefficient of T p : 
16ttS 



X = — : — 7 • (53) 
3Tp{n + a) 

The second term in the flux, with coefficient £ = T%, has been considered to be a purely 
relativistic effect ( Mihalas 1983 ), but there is a classical counterpart when the particles have 
long mean free paths flUhen et al. 200U| ) . We shall see in the next section how these transport 
coefficients, as well as the bulk viscosity, can be deduced from the entropy generation rate. 



4 ENTROPY GENERATION AND BULK VISCOSITY 

To clarify the roles of the various terms in the radiative stress tensor, we compute their 
contributions to the entropy generation. In leading order, this is related to the deviation of 
the stress tensor from its equilibrium value, which we approximate to first order in e as 

AT M1/ = - Tj£ = eTfi , (54) 

where 

Tft = E {1) uV + Ffou" + Ffavf + Pfi . (55) 
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The rate of entropy production is the four-divergence of the entropy flux, S M , whose 
expression as given by Weinberg ([Weinberg 1972) and recalled in Appendix A is 

T££ = AT P [u^ - u M (logT) !V ] . (56) 

For the components of we introduce fl^5|), ( f4"T|) and ( |48"D and obtain, after some simple 
rearrangement s , 



^(log7V + ^ 



Since 5* = aT 4 , we can rewrite 



as 



167T.S 



</(logT), p + -0 



(57) 



(58) 



which we use in the first term of ([57]). To simplify the second term in (|57|), we note that, 
on account of the form of we encounter the expression Up iV T pvpo r u P)Cr . We may use the 
identity h^ vpa = to show that 



u„,,t 



pupa 



Then we find that 

fii/pa 



[U, 



u 



v,p 



U 



p,(7 ■ 



1 2 

^P,<t ~^\Up,u ^v,p ~^tlpv^- 



(59) 



(60) 



yp 



All this adds up to 
T 



where 



16vrS 



3Tp(n + cr) 

32ttS 
3(10k + 9a) 



and 



c 



16ttS 



1 u a T c n 

- H 

3 TO 



(61) 

(62) 
(63) 

(64) 



To see the meaning of these results, we note that for any four- vector v^ 1 , whose compo- 
nents in the comoving frame are (v°,v), we have 

hTv^ = d^ViVj = - || v || 2 < 0. (65) 

Similarly, for any rank-2 tensor w pu , we can show that h pp h ya w IMU w pr j > 0. So we find that 
the rate of entropy generation is never negative as long as all these coefficients are positive. 
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The first term in (|6lD is the entropy generation produced by the drag of the radiative flux, 
which is an effect of radiative thermal conductivity, with due allowance for the acceleration of 
the matter. The second term represents the entropy generation produced by volume changes 
and the coefficient in this term, (, is the radiative bulk viscosity. Finally, the third term is a 
rate of entropy increase caused by shear, with the coefficient of radiative shear viscosity, fi. 

We see from ( |B1D that the bulk viscosity depends on u p T p = T. In depending on a rate, 

our bulk viscosity differs qualitatively from that derived by Weinberg flWeinberg 1972| ) who 

found a bulk viscosity of the form 

I _ fdP\ 
3 ~ Vde~) 



16%S 

Qw 



pnT 

For a radiation-dominated situation, in which (fjf) ~ §, we get (w ~ 0. 



(66) 



The principal reason for the difference between the results (p4|) and (f^) is that Wein- 
berg followed the standard practice in transport theory of using a lower order condition (here 
-^(o) v = 0) t° eliminate the time derivative from the pressure tensor ([Weinberg 1972| ; Pern- 
[stein 1988|) . Thus he obtained S = —AS (j^j & (see Appendix B) and was led to formula 
(|66|) for the bulk viscosity. While this does simplify the derivation of (|66|), it also shrinks its 
domain of validity. In avoiding the use of this condition in our derivation, we obtain a result 
that is capable of describing relatively rapid processes and narrow structures. Moreover, we 
find that the bulk viscous pressure can be large even for the radiation-dominated case. 

It is interesting that with both approaches, as the absorption coefficient gets smaller, the 
bulk viscosity grows larger. This is not so surprising since, for large absorption coefficient, 
the temperature difference between radiation and matter is small. Still, the tendency toward 
large bulk viscosity at small absorption coefficient is noteworthy and is reminiscent of the 
sharply different behaviors of fluids with zero viscosity and those with viscosity tending to 
zero. Of course, since the effect of Compton scattering will be like that of absorption, the 
limit of zero absorption strictly does not apply to the general case. Nevertheless, we may 
expect to find effective heating of expanding media in some situations. 



5 CONCLUSION 

From its very beginnings, the subject of radiative viscosity has been vexatious, as several 
of the founders of modern radiative transfer theory — Milne, Eddington, Jeans — found 
in their attempts to compute the shear viscosity ( Mihalas 1983j ) . It was not until Thomas 
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( [Thomas 1930| ) did the problem relativistically that agreement was reached on the value 
of the shear viscosity, to leading order. What is interesting is that the agreement seems 
to have been reached on purely theoretical grounds since there do not appear to be direct 
measurements confirming Thomas's results. 

In the matter of bulk viscosity, the situation is even more troubling in that the calcu- 
lations are subtler and there are by now many contributors to the growing literature on 
this issue. It is conceivable that disagreements over bulk viscosity may also be settled by 
common consent, but we feel this is unlikely. As we have noted, in previous work, the sim- 
plification of introducing equilibrium approximations into the lower order theory is used, 
in the spirit of the Chapman-Enskog procedure ( |Kogan 1969| ; |Uercignani 1988| ). For gases 
of classical particles, this approximation leads to the conclusion that a simple classical gas 
does not have a bulk viscosity. However, when this restrictive approximation is not made, 
we have found (|Chen et al. 200d[ ), on using the relaxation (or BGKW) model ( [Kogan 1969| ; 



Cercignani 1988|) , that a classical rarefied gas may have an effective bulk viscosity and that 



the gas obeys macroscopic equations that differ from the Navier-Stokes equations when the 
particle mean free paths become long. In the classical case, one has recourse to experiments 
to test the conclusions and these have lent to support to our conclusion. In the radiative 
problem, we may hope for a rough empirical check in the photon density of the present 
universe. This is a calculable quantity, though it is model dependent, and we shall report 
elsewhere on the results obtained for this case (or see ( |(Jhcn 2000 )). 

But this is not the whole story. There has been some previous disagreement on the 
correct form of transport coefficients ( |Andcrson fc Kox 1977| ), especially, the bulk viscosity. 
The calculations are usually done in the frame defined by Eckart and the coefficients are 
identified by comparison with the general form of the stress tensor. But there are usually 
two temperatures in the expression of the stress tensor (for matter and radiation) and 
the difference between them is of the same order as the bulk viscosity. Hence the derived 
bulk viscosity is dependent on which temperature is actually adopted in the stress tensor. 
But with our present method, we needed to go only to first order to get good results, so 
the temperature difference does not figure significantly in the outcome. A calculation of 
the entropy generation rate then leads to an unambiguous determination of the transport 
coefficients. 

Another way to attack this problem is to use the moment method, which was first 
formally developed in kinetic theory by Grad ( |Grad 1963| ); in transfer theory, it may be 
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traced back to Krook's ([Krook 1955|) formulation of Eddington's methods. In modern times, 
the moment method has been elaborated for use in radiative fluid dynamics by several 
authors in attempting to deal with viscous effects, especially nonlocal ones ( |Anderson fc 



Spiegel 1972| ; llhorne 1981| ; |Kato et al. 1993] ; |Uday fc Israel 1982| ; |Struchtrup 1997|) . However, 



in the context of classical kinetic theory, the introduction of higher moments beyond those 
used by Grad does not lead to rapid improvement of the macroscopic description of rarefied 
gases ( phen et al. 20001) . The reason is that the microscopic theory admits only a few slow 
variables, corresponding to the lowest moments of the distribution function. 

Higher moments are fast variables and, as we know from the experience of dynamical 
systems theory, their introduction does not lead quickly to improved representation of the 
dynamics. Rather, it is best to attempt to improve the description in terms of the slow 
variables. Therefore, we have pursued here an approach in which we compute only approxi- 
mations for the energy density, flux and pressure tensor that may be derived by an expansion 
in the style of Thomas. The next step in this approach would be to seek a closure relation 
for the pressure tensor in terms of the radiative variables rather than in terms of a mate- 
rial quantity like the source function S. We shall take this up in a later paper where we 
shall see that, to improve our description in terms of the slow variables, we should include 
their derivatives in the closure relation, since the derivatives of slow variables are also slow 
variables. 

APPENDIX A: DERIVATION OF ENTROPY GENERATION RATE 

Consider a medium made up of N interacting species in a state of near thermal equilibrium. 
Let us represent the one-particle distribution of the zth species as 

y = % + % (ai) 

where 

% = exp[(-f3^ - k) - %] (A2) 

is the equilibrium distribution, with % e = 1,-1 for a Bose-Einstein or a Fermi-Dirac gas, 
respectively. Also, 

h = ^ (A3) 

and h is the chemical potential for the zth species, where % — 1, 2, • • • , N. The corresponding 
4-current and stress tensor are defined as 
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^ = ~ l 9 I -ITP^l l / ln '/ - Ml + '/Ml + * «/)]. (A6) 



Radiative Bulk Viscosity 15 

Yi" = & y ^ , Vr"" =*gj^*f d A p (A4) 

while the corrections are, similarly, 

p - j^p , *rf = l 9] P -jr l f^ 3 P ( A5 ) 

where *g is a normalization constant depending on the nature of the statistics of the particles. 
The entropy four-vector is 
d 3 p 

If we put equation ( |Alj ) into this definition and calculate the entropy generation rate, accu- 
rate to second order in l f\ , we have 

= -h Vi^ - ^ + \ Vif + J T7 . (A7) 

So the total entropy generation rate of the mixture is 

= ^ - E & * T ,r + E M ^ + lJi r] • (as) 

i i i 

Chemical equilibrium leads to J2i %0L ^% = and total energy conservation results in 
J2i ftp l T£ u = 0. With these two constraints, equation (|A8|) reduces to: 

s :; = E H K k + ^ri • (A9) 



With formula (|A9|) we can evaluate the entropy generation rate for radiating media by 
considering only two species in our mixture. The material medium, is one of these and is 
assumed to be in thermal equilibrium in the absence of photons. For this component, f\ = 0. 
The other constituent is the radiation field, whose chemical potential, a, is zero. So we obtain 
for the entropy generation formula of the mixture 

^ = /vr , (aio) 

where T"^ is the correction to the radiative stress tensor. With the definition for /3 V , equation 
( |A10|) readily reduces to the formula reported in the text except that, in this appendix, we 
have not introduced a fiducial e. 



APPENDIX B: THE RELATION BETWEEN S AND 

In order to evaluate the bulk viscosity, we need to express S in terms of 9. This can be done 
by analyzing the zeroth order entropy of the entire system — medium plus radiation — and 
using the conservation laws. We have the thermodynamic relation 

© 2000 RAS, MNRAS 000, 
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Tda = d(ev) + pdv , (Bl) 

where p is the total pressure, a is the entropy per particle, v is the volume per particle and 
e is the total energy density. 

If T and v are chosen to be the independent variables and F is the Helmholtz energy 
per particle, we can write 

dF = -adT - pdv . (B2) 



Since dF is a perfect differential, we have the Maxwell relation 



'da 
dv 



dp N 
dT 



(B3) 



On writing F = ev — Tcr, we find a = k{ev — F). If we put this into the left of the foregoing 



equation and notice that from ( p2|) , 



T 
dF 

dv 



— p, we get 



de 



dp 



v hr +e + p = T\ 



dv 



dT 



We also have the continuity equation 



(B4) 



(B5) 



where h = u^n ^. In terms of v (= the continuity equation can be written as 

v = v?v# = vvPp . (B6) 

If we assume that the bare matter stress tensor satisfies O^j v = with 6^ = eu^u u — 
pW, we find 

e = — (e + p)v/v , (B7) 

where we have used (p6|). Now we are in a position to calculate T in terms of 0(= w^J, 
starting with the identity, 



de 



de \ 



e(v,T) tlt = [ — ) v^ + [ — ) T M . 



dry 



We have 

i= (i) 

Hence 

*-(£. 



e — 



'de\ 
dT) 

(de 
dv 



T . 



v 

T J 



On replacing e from (|B7|), (§|) T from 



(B8) 
(B9) 
(BIO) 

and v from ( |B1j| ) in the above relation, we get 

(Bll) 
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With S = aT A , we find that, in the comoving frame, equation (|B1 1|) reduces to 




S=-4S[%-) 9 . (B12) 
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